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Abstract 

A positive muon is a spin-1/2 particle. Beams of muons with all their spins polarized can be prepared and 
subsequently implanted in various types of condensed matter. The subsequent precession and relaxation of their 
spins can then be used to investigate a variety of static and dynamic effects in a sample and hence to deduce 
properties concerning magnetism, superconductivity and molecular dynamics. Though strictly a lepton, and behaving 
essentially like a heavy electron, it is convenient to think of a muon as a light proton, and it is often found with 
a captured electron in a hydrogen-like atom known as muonium. This article outlines the principles of various 
experimental techniques which involve implanted muons and describes some recent applications. The use of muons 
in condensed matter physics has shed new light on subjects as diverse as passivation in semiconductors, frustrated 
spin systems, vortex lattice melting, and quantum diffusion of light particles. 

Originally published in Contemporary Physics 40, 175, (1999) 



1 Introduction 

Condensed matter physics is concerned with the prop- 
erties of solids, liquids and various intermediate states 
of matter such as colloids (Chaikin and Lubensky 1995). 
From a fundamental viewpoint, all matter is made up of 
two main types of constituent particles, quarks and lep- 
tons (see Figure |l|). The universe is composed of only 
these particles in addition to the gauge bosons (such as 
photons) which mediate forces between particles. Pro- 
tons and neutrons are each made up of three quarks 
while electrons are thought to be indivisible and are 
members of a class of particle called leptons. Everyday 
matter, made up of protons, neutrons and electrons, 
is the main concern of condensed matter physics and 
therefore the subject arises from particles in the first 
column (or 'generation') of the standard model (Fig- 
ure |l]). However, more exotic particles do exist which 
are either more unusual leptons or combinations of more 
exotic quarks; they are all short-lived and have never 
been found to be relevant in condensed matter physics 
with one single exception: the muon. 

The muon is found in nature; it is the dominant con- 
stituent of the cosmic rays arriving at sea-level and 
a few will be hitting you each minute as you read 
this article. The fact that so many muons make it 
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Figure 1: Quarks (u=up, c?=down, s=strange, 
c=charm, i=top, fo=bottom) and leptons (e=electron, 
/x=muon, r=tau, and their associated neutrinos z/g, v^, 
Vr . Everyday matter is composed of particles from the 
first column in this table. 
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through the atmosphere is a weU-known consequence 
of time-dilation effects in special relativity (Rossi and 
Hall 1941, Bailey et al. 1977). However, at a number 
of locations in the world, intense beams of muons are 
prepared artificially for research in condensed matter 
physics. In this article I will try and explain why this 
is done and what can be learned by firing this rather 
exotic 'second-generation' particle into various types of 
matter. I will describe how implanted muon studies 
have shed new light on problems in condensed mat- 
ter physics, including those concerning semiconductors, 
magnetism, superconductors, and quantum diffusion. 
A key difference between this technique and those in- 
volving neutrons and X-rays is that scattering is not in- 
volved. Neutron diffraction methods use the change in 
energy and/or momentum of a scattered neutron to in- 
fer details about a particular sample. In contrast muons 
are implanted into a sample and reside there for the 
rest of their lives. The muons themselves never emerge 
again. It is the positrons into which they decay that are 
released from the sample and yield information about 
the muons from which they came. The experimenter 
plays the role of a coroner at an inquest; the evidence 
of the decay particles given up at death allows one to 
infer the nature of the muon's life. 

2 The discovery of the muon 

Although the discovery of the muon is usually consid- 
ered as a single event which is credited to Neddermeyer 
and Anderson because of their work in 1936, the story is 
actually more complicated and one might say that the 
muon slowly "emerged" over almost a half-century. Fol- 
lowing the discovery of radioactivity in 1899, it was ob- 
served that electrometers, instruments capable of mea- 
suring the ionization produced by radioactivity, dis- 
charged even when there was no obvious radioactive 
source. A Jesuit priest called Thomas Wulf noticed in 
1910 that this effect was more pronounced at the top 
of the Eiffel tower than at the bottom (Wulf 1910, see 
also Xu and Brown 1987). That the rate of discharge 
of an electrometer was an increasing function of alti- 
tude was demonstrated beyond doubt by Victor Hess 
in an intrepid series of ballooning experiments in the 
period 1911-12 (Hess 1912). He showed that electrom- 
eters discharged between three and five times faster at 
an altitude of 5000 m than they did at sea level. Such 
experiments were later extended to very high altitudes 
by the use of remote-controlled electrometers (in which 
the experimenter did not actually accompany their ex- 
periment in the balloon); this advance was achieved by 
Robert Millikan who was initially sceptical about Hess' 
work but became rapidly convinced of its importance 



(Xu and Brown 1987). Millikan's name for the effect, 
cosmic rays, gained universal acceptance, but it was 
Hess who was awarded the Nobel prize for the discov- 
ery in 1936. The major constituent of cosmic rays at 
ground level turned out to be the muon, but the true 
identity of the muon remained hidden for some time. 

Millikan's group at CalTech began to perform exper- 
iments on cosmic rays at ground level by bending the 
rays in magnetic and electric fields. This work led to 
the discovery of the positron by Carl Anderson in 1932 
(Anderson shared the Nobel prize in 1936 with Hess) 
and then the muon (initially called the mesotron, and 
later the mu-meson) was identified in cosmic rays by 
Seth Neddermeyer and Carl Anderson in 1936. They 
measured the mass of the muon, which turned out to 
be roughly one-ninth of the proton mass, and agreed 
extremely well with the predictions of the Japanese 
physicist Hideo Yukawa for a particle to mediate the 
strong force and hold the nucleons in the nucleus to- 
gether. This result can be understood by using Heisen- 
berg's uncertainty principle AiJAt ~ h for this virtual 
particle. Ai? is the the rest mass energy of Yukawa's 
particle, the energy which has to be "borrowed", and 
Af is the time for it to cross the nucleus, the time for 
which it has to be "borrowed" ; by assuming a speed 
~ c for the particle one arrives at a mass ~100 MeV, 
the mass found for the muon. The muon was therefore 
assumed to be the particle which holds the nucleus to- 
gether. It took some technically difficult experiments 
to show that this assumption was in error and that the 
startling agreement with Yukawa's predictions was mis- 
leading. Marcello Conversi, Ettore Pancini and Oreste 
Piccioni, working in a basement in Rome during and 
after the Second World War, tried to implant cosmic 
ray muons in matter and then measure their lifetime 
(Conversi et al. 1945, 1947). They found that positive 
muons implanted in anything always live on average for 
2.2 /xs (the lifetime in vacuum); the lifetime of negative 
muons on the other hand depends on the atomic num- 
ber Z of the material into which it is implanted; they 
measured 2.2 /us for carbon but found 0.07 /xs for lead. 
They reasoned that if muons really mediated the strong 
force they should be gobbled up much more quickly in 
all materials. It appeared that muons interacted only 
rather weakly with matter and the only effect that was 
observed occurred for negative muons and turned out 
to be /x~ capture: 

A negative muon is attracted by atomic nuclei and since 
its mass is much greater than that of an electron, it 
readily displaces an electron from an atom and rapidly 
drops down to the Is state. Prom there it either de- 
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1:1: 2 : 3 

Figure 2: The muon can be considered to be the light- 
est member of the family of particles which include the 
proton, deuteron and triton. The approximate ratio 
of their masses are as shown. The muon thus extends 
the mass range of isotope substitution available to the 
experimenter. 



cays to an electron or undergoes capture as above. The 
propensity to undergo capture depends on the atomic 
number Z (approximately as the fourth power of Z) 
through the Bohr radius of the /i~ orbit. 

Yukawa's particle turned out to be the pion, discov- 
ered in cosmic rays by Powell's group in 1947 (Lattes et 
al. 1947), which decays into a muon after only 0.026 ^s, 
making it harder to spot in cosmic rays. Thus the muon 
was not a meson (a term now reserved for quark-anti- 
quark pairs like the pion) but a lepton, a heavy electron 
with no internal structure. 

For experiments in condensed matter physics it is 
mainly the positive muon which is used. The negative 
muon fi~ , which implants close to an atomic nucleus, is 
generally much less sensitive to the phenomena of inter- 
est to condensed matter physicists (magnetism, super- 
conductivity, etc) than the site of the implanted positive 
muon which sits well away from nuclei in regions of 
large electron density. Condensed matter physics is es- 
sentially the physics of electrons, rather than nuclei, so 
that the best place you can put your test particle is in 
the electron cloud. In fact even though the muon is a 
lepton (see Figure |]) and therefore essentially a heavy 
electron, for our purpose it is more useful to consider it 
as a light proton (see Figure H) . 

Some properties of the electron, muon and proton are 
tabulated in Figure |. The mass of the muon is inter- 
mediate between that of the electron and the proton, 
and thus so are its magnetic moment and gyromagnetic 
ratio. The latter is the constant of proportionality be- 
tween angular momentum and magnetic moment. 
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Figure 3: Properties of the electron, muon and proton. 

3 Muon production, decay and 
implantation 

Cosmic rays provide a major source of muons [roughly 
one muon arrives vertically on each square centimetre 
of the earth's surface every minute (Caso et al. 1998)] 
but, since the 1950s, experiments on muons have needed 
higher intensities and have required the development 
of accelerators. High energy proton beams (produced 
using synchrotrons or cyclotrons) are fired into a target 
(usually graphite) to produce pions via 

p + p ^ + p + n, 

and the pions decay into muons: 

where Vfj^ is a muon-neutrino. The pion decay is a two- 
body decay and is therefore particularly simple. For 
example, consider the pions which are produced at rest 
in the laboratory frame. To conserve momentum, the 
muon and the neutrino must have equal and opposite 
momentum. The pion has zero spin so the muon spin 
must be opposite to the neutrino spin. One useful prop- 
erty of the neutrino is that its spin is aligned antiparallel 
with its momentum (it has negative helicity), and this 
implies that the muon-spin is similarly aligned. Thus by 
selecting pions which stop in the target (and which are 
therefore at rest when they decay) one has a means of 
producing a beam of 100% spin-polarized muons. This 
is the method most commonly used for producing muon 
beams for condensed matter physics research, though 
other configurations are in use (Brewer 1994). 

The muons are stopped in the specimen of interest 
and decay after a time t with probability proportional 
to e~*/'^f' where = 2.2 /xs is the lifetime of the muon. 
The muon decay is a three body process 

^+ e+ + Z^e + V^i 

and so the energy of the positron e~^ (which is the only 
particle produced in this reaction that we have a sensi- 
ble hope of reliably detecting) may vary depending on 



Spin-polarized muons in condensed matter physics 



S. J. Blundell 



4 



how momentum is distributed between the three par- 
ticles (subject to the constraint that the total vector 
momentum will sum to zero, the initial momentum of 
the stopped muon). The decay involves the weak in- 
teraction and thus has the unusual feature of not con- 
serving parity (Garwin et al. 1957). This phenomenon 
(which also lies behind the negative helicity of the neu- 
trino) leads to a propensity for the emitted positron to 
emerge predominantly along the direction of the muon- 
spin when it decayed. This can be understood by con- 
sidering the mirror image of the muon decay process 
(Figure |). 

The angular distribution of emitted positrons is 
shown in Figure ^ for the case of the most energeti- 
cally emitted positrons. In fact positrons over a range 
of energies are emitted so that the net effect is some- 
thing not quite as pronounced, but the effect neverthe- 
less allows one to follow the polarization of an ensemble 
of precessing muons with arbitrary accuracy, providing 
one is willing to take data for long enough. 

Muons are implanted into the sample with an energy 
which is at least 4 MeV. They lose energy very quickly 
(in 0.1-1 ns) to a few keV by ionization of atoms and 
scattering with electrons. Then the muon begins to un- 
dergo a series of successive electron capture and loss 
reactions which reduce the energy to a few hundred eV 
in about a picosecond. If muonium is ultimately formed 
then electron capture ultimately wins and the last few 
eV are shed by inelastic collisions between the muo- 
nium atom and the host atoms. All of these effects are 
very fast so that the muon (or muonium) is thermalized 
very rapidly. Moreover the effects are all Coulombic in 
origin and do not interact with the muon-spin so that 
the muon is thermalized in matter without appreciable 
depolarization. This is a crucial feature for muon-spin 
rotation experiments. One may be concerned that the 
muon may only measure a region of sample which has 
been subjected to radiation damage by the energetic 
incoming muon. This does not appear to be a prob- 
lem since there is a threshold energy for vacancy pro- 
duction, which means that only the initial part of the 
muon path suffers much damage. Beyond this point of 
damage the muon still has sufficient energy to propa- 
gate through the sample a further distance thought to 
be about 1 /im, leaving it well away from any induced 
vacancies (Chappert 1984). 




Figure 4: Top: molecules in a box collide with the 
walls. These collisions do not violate parity so that 
both the process shown on the left, and its mirror im- 
age on the right, could be observed in nature. Bottom: 
the same is not true for the process of muon decay. The 
direction of the muon-spin is reversed in the mirror so 
that the positrons are emitted predominantly in a direc- 
tion opposite to that of the muon-spin. The violation 
of parity means that in our universe only the process 
on the left-hand side of the diagram is ever observed. 



4 Spin precession and relaxation 

In a magnetic field B the muon-spin precesses with an- 
gular frequency oj^ given by uj^ = j^B where 7^ — 
ge/2m^ is the gyromagnetic ratio for the muon. This 
is known as Larmor precession. The field-dependent 
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Figure 5: The angular distribution of emitted positrons 
with respect to the initial muon-spin direction. The fig- 
ure shows the expected distribution for the most ener- 
getically emitted positrons. 



precession frequencies for the muon, electron and pro- 
ton are shown in Figure ^. The highest frequencies are 
associated with the lightest particle, the electron, and 
the lowest with the proton. For usual laboratory mag- 
netic fields this explains why ESR (electron-spin reso- 
nance) is typically performed at microwave frequencies 
while NMR (nuclear magnetic resonance) uses radio- 
frequencies. In both of these techniques resonance oc- 
curs when the precession frequency matches the reso- 
nance frequency. Muon-spin rotation (/iSR) is associ- 
ated with frequencies intermediate between NMR and 
ESR but unlike those resonance techniques, no electro- 
magnetic field is necessary since the precessing muon 
can be followed directly. [Muon-spin resonance experi- 
ments can however be performed (Brewer 1994) but a 
discussion is outside the scope of this article.] 

A schematic diagram of the experiment is shown in 
Figure 0(a). A muon, with its polarization aligned an- 
tiparallel to its momentum, is implanted in a sample. 
(It is antiparallel because of the way that it was formed, 
see above, so the muon enters the sample with its spin 
pointing along the direction from which it came.) If the 
muon is unlucky enough to decay immediately, then it 
will not have time to precess and a positron will be 
emitted preferentially into the backward detector. If it 
lives a little longer it will have time to precess so that 
if it lives for half a revolution the resultant positron 
will be preferentially emitted into the forward detector. 
Thus the positron beam from an ensemble of precess- 
ing muons can be likened to the beam of light from a 
lighthouse. 

The time evolution of the number of positrons de- 
tected in the forward and backward detector is de- 
scribed by the functions Npit) and Nsit) respectively 
and these are shown in Figure 0(b). Because the muon 
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Figure 6: The Larmor precession frequency / in MHz 
(and the corresponding period r = 1//) for the elec- 
tron, muon and proton as a function of applied mag- 
netic field B. 
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Figure 7: Schematic illustration of a /iSR experiment, 
(a) A spin-polarized beam of muons is implanted in a 
sample S. Following decay, positrons are detected in 
either a forward detector or a backward detector B. 
If a transverse magnetic field H is applied to the sample 
as shown then the muons will precess. (b) The num- 
ber of positrons detected in the forward and backward 
detectors, (c) The asymmetry function. 
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Figure 8: Zero- field muon spin rotation frequency in the 
organic ferromagnet p-NPNN (Blundell et al. 1995). 

decay is a radioactive process these two terms sum to an 
exponential decay. Thus the time evolution of the muon 
polarization can be obtained by examining the normal- 
ized difference of these two functions via the asymmetry 
function A{t), given by 



A[t) 



NB{t)~NF{t) 

NB{t) + NF{ty 



(1) 



and is shown in Figure ^(c). 

This experimentally obtained asymmetry function 
has a calculable maximum value, Amax, for a partic- 
ular experimental configuration which depends on the 
initial beam polarization (usually very close to 1), the 
intrinsic asymmetry of the weak decay, and the effi- 
ciency of the detectors for positrons of different ener- 
gies, and usually turns out to be around Amax ~0.25. 
The function can be normalized to 1, in which case 
it expresses the spin autocorrelation function of the 
muon, G{t) — A(t)/Amax, which represents the time- 
dependent spin polarization of the muon. 

A magnetic field does not need to be applied for 
the muons to precess if the sample has its own mag- 
netic field. Figure ^ shows muon-spin rotation data for 
muons implanted into an organic ferromagnet (Blundell 
et al. 1995). This material magnetically orders only at 
very low temperatures (Tc = 0.67 K) so that exper- 
iments must be carried out in a dilution refrigerator. 
This presents no problems for the muon which passes 
through the windows of the cryostat and implants in 
the sample, revealing the internal magnetisation. As 
the sample is warmed, the frequency of oscillations de- 
creases as the internal field decreases until it is above 
the Curie temperature and no oscillations can be ob- 
served, only a weak spin relaxation arising from spin 
fiuctuations. 



This is an example of an experiment with no applied 
magnetic field (in fact a small magnetic field was ap- 
plied to compensate the effect of the earth's field) . Very 
often magnetic fields are applied to the sample either 
perpendicular or parallel to the initial muon-spin di- 
rection. The perpendicular (or transverse) case causes 
the muon to precess in the applied magnetic field and 
any dephasing in the observed oscillations is evidence 
for either an inhomogeneous internal field distribution 
or spin-spin (T2 in the language of NMR) relaxation. 
The parallel (or longitudinal) case will be described in 
more detail later and does not lead to spin precession, 
but spin relaxation. This can be due to inhomogeneous 
field distributions or spin-lattice (Ti) relaxation pro- 
cesses. 

Muon experiments can be performed in two differ- 
ent ways depending on the time structure of the muon 
beam. If the muon beam is continuous (CW, or con- 
tinuous wave see Figure §(a)), then muons arrive at 
the sample intermittently. When the muon enters the 
experiment it must itself be detected to start a clock. 
When the positron is detected in either the forward or 
backward detectors, the clock is stopped. If a second 
muon arrives before the first one has decayed then one 
has no way of knowing whether a subsequently emitted 
positron came from the first or second muon, so this 
event must be disregarded. Sophisticated high-speed 
electronics and a low incident muon arrival rate are 
needed. Alternatively one can use an electrostatic de- 
flector triggered by the detectors to ensure no muons 
enter the experiment until the current implanted muon 
decays. 

These complications are circumvented with a pulsed 
muon beam (Figure gb)). In this case a large number 
of muons arrive in a very intense pulse so there is no 
need to detect when each muon arrives. The detection 
of positrons is then made and each event is timed with 
respect to the arrival of the pulse. A typical dataset 
contains several million detected positrons so that an 
appreciable number of muons (the fraction is given by 
g-20/2.2 ^0.01%) live for 20 /zs or longer. Long-lived 
muons are difficult to measure with CW beams; the 
the arrival of the next muon tends to interrupt the first 
muon which has outstayed its welcome! Nevertheless 
the long-lived muons can be accurately detected at a 
pulsed source. Unfortunately this method also suffers 
from a drawback which is that the muon pulse has a 
finite width, r^,, which results in a slight ambiguity in 
all of the timing measurements and leads to an upper 
limit on precession frequencies which can be measured. 
By the standard of other unstable elementary particles, 
the muon is comparatively long-lived with a lifetime of 
T =2.2 /IS. CW muon beams are operated at the Paul 
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Figure 9: Schematic illustration of the two types of 
muon beam (a) continuous wave (CW) and (b) pulsed. 
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Figure 10: The Breit-Rabi diagram for isotropic muo- 
nium with a hyperfine constant A ~ 4463 MHz appro- 
priate for free (vacuum) muonium. Also shown are the 
splittings between pairs of levels over a larger range of 
magnetic field. Levels 1 and 2 cross at a field of ~ 16.5 T 
for free muonium. 

Scherrer Institute in Switzerland and at TRIUMF in 
Canada. Pulsed muon beams are used at KEK in Japan 
and at ISIS (the spallation source at the Rutherford 
Appleton Laboratory) in the UK. 

5 Muonium 

Depending on its chemical environment, the muon can 
thermalize and pick up an electron and form a neutral 
atomic state called muonium (abbreviated Mu= /i+e^) 
which is an analogue of atomic hydrogen. In muonium 
the electronic spin and the muon-spin are coupled by 
a hyperfine interaction which we will initially assume 
is isotropic. This leads to two energy levels, a lower 
triplet state and a higher singlet state. In a magnetic 
field the triplet levels split and the energy levels move 
as shown in the Breit-Rabi diagram in Figure |l^. 
Isolated hydrogen or muonium is one of the simplest 



defects in a semiconductor. The electronic structures 
of analogous hydrogen and muonium centres are ex- 
pected to be identical, apart from vibrational effects 
reflecting their difference in mass. If muons are im- 
planted into a semiconductor like silicon, neutral muo- 
nium tends to form at the tetrahedral sites (labelled T 
in Figure |ll|) and can diffuse rapidly between such sites. 
The electronic state is isotropic and in a transverse-field 
experiment the transitions between energy levels can 
be observed as precession frequencies Vij equal to the 
splitting between energy levels [Ei — Ej)/h (as shown 
in Figure |l0|). The strength of the hyperfine interac- 
tion in semiconductors is usually close to half that of 
the vacuum value (the reduction is due to some ad- 
mixture of the electron spin density with surrounding 
atoms). However, a substantial fraction of neutral muo- 
nium is also found in a most unexpected place, wedged 
into the centre of a stretched Si-Si bond (the site is 
labelled BC in Figure |ll| for 'bond-centre', see Patter- 
son 1988 for a review of muonium states in semicon- 
ductors). This state is extremely immobile, and sur- 
prisingly turns out to be the thermodynamically more 
stable site. Its hyperfine coupling is much lower than 
that of the tetrahedral state, typically less than 10% 
of the vacuum value. Furthermore the coupling is very 
anisotropic, with axial symmetry about the (111) crys- 
tal axis (i.e. along the Si-Si bonds), so that the energy 
levels behave in a different manner to that indicated 
in Figure |l0[ These states have rather interesting dy- 
namics and can undergo charge and spin exchange pro- 
cesses, cycling rapidly between positive and negative 
charge states via interaction with conduction electrons 
(Chow et al. 1994). The characterisation of all these 
states is important because it is found that atomic hy- 
drogen is present in most semiconductors and is able to 
passivate (i.e. deactivate) the dangling bonds in amor- 
phous silicon allowing it to show semiconducting prop- 
erties. Hydrogen is inevitably present in all semicon- 
ductors, often becoming incorporated during material 
production from hydride gases or during etching, but 
the low concentration makes direct spectroscopic stud- 
ies very difficult. Using muonium (albeit in the ultra- 
dilute limit!) as an analogue for hydrogen has therefore 
been a promising method of obtaining a great deal of 
spectroscopic information concerning this problem. 

Muonium states can also be formed in many chemical 
systems and allow a unique form of radical spectroscopy 
(Roduner 1993). Muonium adds to unsaturated bonds 
to form muonated free radicals. For example, addi- 
tion to benzene (CgHg) leads to the muonated cyclo- 
hexadienyl radical (CeHgMu). The advantage here is 
that one can work with concentrations down to just one 
muonated radical at a time in an entire sample. In con- 
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Figure 11: The crystal structure of silicon, showing the 
possible muonium sites (T=tetrahedral site, BC=bond- 
centre site). 
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Figure 12: Repolarization of the initial muon polar- 
ization. 

trast ESR detection needs ~ 10^^ radicals in a cavity, 
forbidding measurements at high temperatures where 
the radicals become mobile and terminate by combi- 
nation. The technique has been applied to radicals in 
various environments (Roduner 1993) including those 
absorbed on surfaces (Reid et al. 1990). 

In all cases the muonium can be studied by mea- 
suring precession signals in an applied magnetic field, 
or by using a technique known as repolarization (Fig- 
ure |lj). In this latter method a longitudinal magnetic 
field is applied to the sample, along the initial muon- 
spin direction, and as the strength of the magnetic field 
increases, the muon and electron spins are progressively 
decoupled from the hyperfine field. For isotropic muo- 
nium, half of the initial polarization of implanted muons 
is lost because of the hyperfine coupling, but this is re- 
covered in a sufficiently large applied field (Figure 12), 
allowing an estimate of the strength of the hyperfine 
field. If the hyperfine coupling is anisotropic then this 
affects the form of the repolarization curve, allowing 



further information to be extracted (Pratt 1997). 

At very high fields, when the electron and muon spins 
are completely decoupled, the initial muon polarization 
is preserved. However, at certain values of the magnetic 
field a level-crossing of energy levels may occur. In fact, 
there is one at very high field in vacuum muonium (see 
Figure [Tc|(b)). Interactions between the muon-clcctron 
system and nuclei in the host material cause the pure 
Zeeman states to mix near these level-crossings, thereby 
avoiding the crossing, and this can lead to a loss of po- 
larization. This provides a further technique to extract 
hyperfine coupling strengths and also measurements of 
quadrupole couplings (Cox 1987, Roduner 1993, Brewer 
1994). It turns out that these avoided- level-crossing res- 
onances are also a sensitive probe of dynamics and this 
technique has been applied to a number of studies of 
molecular motion (Roduner 1993). 

An interesting situation occurs in the case of Cgo in 
which muonium can implant inside the buckyball cage 
(this state is called endohedral muonium). The un- 
paired electron part of the muonium greatly enhances 
the sensitivity to scattering from conduction electrons 
so that this state is extremely useful for studying alkali- 
fuUeride superconductors (MacFarlane et al. 1998). It 
is also possible to form a muonium radical by external 
addition, essentially muonium attacking the outside of 
a buckyball, breaking a double bond and ending up 
covalently bonded to a single saturated carbon atom. 
This centre is very sensitive to the molecular dynamics 
of the local environment and has been used to extract 
the correlation time for molecular reorientation (Kiefl 
et al. 1992). 

In metallic samples the muon's positive charge is 
screened by conduction electrons which form a cloud 
around the muon, of size given by a Bohr radius. Thus 
rather than muonium, is the appropriate particle to 
consider in a metal. (The endohedral muonium found 
in alkali fuUeride superconductors is the only known 
example of a muonium state in a metal.) In insulators 
and semiconductors screening cannot take place so that 
the muon is often observed in these systems either as 
muonium or is found to be chemically bound to one of 
the constituents, particularly to oxygen if it is present. 
Isotropic muonium states are found in many semicon- 
ducting and insulating systems. The value of the hy- 
perfine coupling strength is close to that for vacuum 
(free) muonium if the band gap is large. For materials 
with smaller band gap the hyperfine coupling is lower 
reflecting the greater delocalization of the electron spin 
density on to neighbouring atoms (Cox 1987). 
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6 Muons and magnetism 

Muons are ideally suited to studying problems in mag- 
netism. Implanted muons in magnetically ordered ma- 
terials precess in the internal magnetic field and directly 
yield signals proportional to that magnetic field. In 
this respect the muon behaves as a microscopic mag- 
netometer. The very large magnetic moment of the 
muon makes it very sensitive to extremely small mag- 
netic fields (down to ~ 10^^ T) and thus is very useful 
in studying small moment magnetism. It is also valu- 
able in studying materials where the magnetic order is 
random or of very short range. Since muons stop uni- 
formly throughout a sample, each signal appears in the 
experimental spectrum with a strength proportional to 
its volume fraction, and thus the technique is helpful in 
cases where samples may be multiphase or incompletely 
ordered. Because no spatial information is obtained (in 
contrast to diffraction techniques) single crystal sam- 
ples are not essential (though they can be useful in 
certain cases) and experiments can often provide in- 
formation on the magnetic order of certain materials 
where conventional magnetic neutron diffraction can- 
not be simply performed. 

The most straightforward application of the tech- 
nique is to ferromagnets and antiferromagnets in which 
cases the muon is used to follow the temperature de- 
pendence of the internal field and to extract critical 
exponents (Schenck and Gygax 1995). An example of 
this type of experiment has already been presented in 
Figure |[ The internal field at the muon site is due 
to a sum of the applied field (if used), the dipolar 
and demagnetization fields (which can be calculated 
from the magnetization) and the hyperfine field induced 
by the applied field. To extract quantitative informa- 
tion from /iSR experiments it is necessary to know the 
muon-site and this can in some cases hinder the search 
for a straightforward interpretation of the data. Usu- 
ally there are a small set of possible interstitial sites 
which the muon can occupy and in favourable circum- 
stances only one will be consistent with the observed 
data. Nevertheless the technique has been widely ap- 
plied to magnetic materials (Schenck and Gygax 1995, 
Dalmas de Reotier and Yaouanc 1997) and has found 
great applicability to the study of heavy fermion sys- 
tems (Amato 1997). These latter materials are based on 
rare-earth and actinide elements and show strong elec- 
tronic correlations between localized / moments and 
conduction electrons. They exhibit a subtle compe- 
tition between the Kondo effect (by which magnetic 
moments are 'mopped up' by the screening spin po- 
larization of conduction electrons) and the magnetic 
Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions 




Figure 13: Muon-spin precession with a magnetic field 
B applied at an angle 6. 

(which couple magnetic moments through the conduc- 
tion electrons). Depending on the relative strength of 
these effects, magnetic order can be observed at low 
temperature but it is often incommensurate or random 
and can be associated with very small static moments. 
The sensitivity of the muon to tiny fields and the ab- 
sence of a muon quadrupolar moment (which can com- 
plicate analogous NMR experiments) has led to muons 
being extensively utilised in this field (Amato 1997). 

To understand the ability of the muon to study ran- 
domness and dynamics in magnetism it is helpful to 
consider further some aspects of spin precession. If the 
local magnetic field at the muon-site is at an angle of 
6 to the initial muon spin-direction at the moment of 
implantation, the muon-spin will subsequently precess 
around the end of a cone of semi-angle 9 about the mag- 
netic field (Figure [T^). The normalized decay positron 
asymmetry will be given by 



G{t) = cos^ 9 + sin^ 9 cos(7pBt). 



(2) 



If the direction of the local magnetic field is entirely 
random then averaging over all directions would yield 



(3) 



If the strength of the local magnetic field is taken from 
a Gaussian distribution of width centred around 

zero, then a straightforward averaging over this distri- 
bution gives 



Git) 



(4) 



a result which was first obtained by Kubo and Toy- 
abe (Kubo and Toyabe 1967) as an entirely theoreti- 
cal exercise. This relaxation function is illustrated in 



Figure 
Figure 



14 



(b). Its origin is indicated schematically in 
14(a) which shows a number of curves of equa- 
tion ^ for different values of the internal field B. Ini- 
tially they all do roughly the same thing (i.e. fall from 
1 to a minimum value and then increase) but after a 
short time they dephase with respect to each other. 
Hence their average, the Kubo and Toyabe relaxation 
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function, would be expected to fall from unity to a min- 
imum and then recover to an average value, in this case 
to one-third. 

This relaxation function is often observed experimen- 
tally, for example in experiments on copper at 50 K with 
zero applied field (Luke et al. 1991). At this temper- 
ature the muon is stationary and precesses in the field 
due to the neighbouring nuclear dipoles which are ran- 
domly orientated with respect to each other and give 
rise to a field distribution each component of which is 
Gaussian distributed about zero. (At higher temper- 
atures this is not seen because of thermally activated 
muon hopping while at lower temperatures an effect 
known as quantum diffusion may occur, see section ^) 

If the form of the distribution of internal fields was 
different then this would affect the form of the ob- 
served muon-spin time evolution. For example, in a 
material with a spin-density wave which is incommen- 
surate with the crystal lattice, there will be a sinusoidal 
modulation of the internal field which the muons will 
randomly sample. In this case one finds that the muon- 
spin relaxation follows a Bessel function (Major et al. 
1986, Amato 1997). Such an effect has been observed 
in chromium (Major et al. 1986) and also in an organic 
system, (TMTSF)2PF6, which exhibits a spin-density 
wave ground state (Le et al. 1993). 

If there is an almost uniform static internal field in 
the sample, but there is a slight variation from site to 
site, different muons will precess at slightly different 
frequencies and become progressively dephased so that 
the oscillations in the data will be damped. If the field 
varies a great deal the damping could be so large that 
no oscillations can be observed. However this effect 
could also be caused by fluctuations either of the in- 
ternal held, because of some intrinsic property of the 
sample, of the muon's position, or because of muon dif- 
fusion. One method of distinguishing between these 
possibilities is to apply a magnetic field in the longitudi- 
nal direction, parallel to the initial muon-spin direction. 
For example, this modifies the Kubo-Toyabe relaxation 
function as shown in Figure ^(c), causing the "i-tail" 
to increase since in this case the muons precess in both 
the internal field and the applied field. Since the mag- 
netic field Bl is applied along the initial muon-spin 
direction then in the limit of large Bl the muon-spin 
is held constant and does not relax from a value near 
unity. This technique is useful to distinguish the ef- 
fects of inhomogeneous line broadening (a distribution 
of static internal fields as considered above) and fluctu- 
ations because the two have very different behaviours 
in a longitudinal field (Hayano et al. 1979). 

The effects of muon hopping on the relaxation are 
shown in Figure fsl The different traces are for dif- 
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Figure 14: (a) The time evolution of the muon-spin 
polarization for equation ^ with different values of the 
magnitude of the local field \B\. (b) The averaging 
of terms from (a) yields the Kubo-Toyabe relaxation 
function (equation ^ with its characteristic dip and 
recovery to a value of 1/3. (c) The effect of a longitu- 
dinal magnetic field B^ ■ The time is measured in units 
of A^^ and the longitudinal field values shown are in 
units of 
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Figure 15: The relaxation function for a muon hop- 
ping at rate v. After each hop the value of the internal 
field is taken from a Gaussian distribution about zero 
with width The curve for v — corresponds to 

the zero- field Kubo-Toyabe relaxation function of equa- 
tion |4[ 

ferent hopping rates v. When v = we recover the 
zero-field Kubo-Toyabe curve of equation ^. For fast 
hopping the relaxation of the muon-spin becomes dom- 
inated by the hopping process and the relaxation is ex- 
ponential. The relaxation rate goes down as the dy- 
namics get faster. This effect is known as motional 
narrowing because it is essentially identical to the ef- 
fect of the same name in NMR spectroscopy. Motional 
effects narrow NMR linewidths which are measured in 
the frequency domain; consequently in the time domain 
(in which data from /iSR experiments are usually con- 
sidered) this corresponds to a motional widening and a 
reduction in relaxation. 

For slow hopping, very little effect is observed at very 
short times times but a large sensitivity to weak hop- 
ping is observed in the i-tail which is observed at longer 
times. This sensitivity to slow dynamics via the be- 
haviour of the tail of the relaxation function observed 
at long times permits a measurement of dynamics over 
a very large range of time scale. 

If a longitudinal magnetic field is applied then it has 
a large effect on the muon relaxation if the dynamics 
are weak but much less of an effect if the dynamics are 
fast. Thus by a careful combination of zero-field and 
longitudinal-field experiments the nature of the internal 
field distribution can be extracted. 

Such an analysis has been done for many systems 
in which the issue of the presence of dynamics and/or 
magnetic order has profound consequences. A good ex- 
ample is that of the spin ladder cuprates (Kojima et 
al. 1995). Sr2Cu406 has a crystal structure such that 
Cu^+ spins are arranged in weakly coupled 2-leg lad- 



ders. /iSR experiments show no magnetic order down 
to very low temperatures. Sr4Cu60io is a 3-leg lad- 
der and /LtSR finds static magnetic order appearing at 
^^52 K. This is consistent with theoretical predictions 
of a non-magnetic ground state in even-leg-number sys- 
tems (due to the formation of spin-singlets between 
pairs of spins on each rung) but a magnetic ground 
state in odd-leg-number systems (Kojima et al. 1995). 
Similar types of experiments have been performed on 
SrCr8Ga40i9, a material which has magnetic moments 
lying in an arrangement known as a Kagome lattice 
which is highly frustrated in the sense that there is no 
unique ground state. These measurements have demon- 
strated the presence of dynamic spin fiuctuations which 
persist down to at least 0.1 K (Uemura et al. 1994). 
These strong dynamics which remain at very low tem- 
peratures are characteristic of a spin-liquid state in 
which a small number of unpaired spins migrate in a sea 
of singlet pairs. Some compounds with one-dimensional 
chains of antiferromagnetically coupled spins can show 
an effect known as the spin-Peierls transition. Below a 
transition temperature the chain dimerises, pairing al- 
ternate spins into singlets and opening up a gap in the 
excitation spectrum, drastically changing the dynam- 
ics. /iSR studies have been performed on both inorganic 
and organic systems which show this effect (Lappas et 
al. 1994, Blundell et al. 1997). 

A celebrated example of a frustrated spin system is a 
spin glass, such as the dilute alloy spin glass prepared 
by dissolving small concentrations of Mn in a Cu ma- 
trix. In these system dilute magnetic impurities couple 
via an RKKY exchange interaction which alternates in 
sign as a function of distance. Because the magnetic 
impurities are present at random, these materials can- 
not show long range order but have built-in frustration. 
When cooled one observes a slowing down of spin fluctu- 
ations and a divergence in the correlation time between 
Mn spin fluctuations at the spin glass temperature. Be- 
low this temperature, a static component of the local 
fleld is observed with muons, corresponding to some de- 
gree of spin freezing, with each Mn spin having its own 
preferred orientation, although fluctuating around this 
(Uemura et al. 1985). Above the spin glass temperature 
muon-spin relaxation measurements have been used to 
follow the spin glass dynamics and to directly extract 
the form of the autocorrelation function of the spins 
(CampbeU et al. 1994, Keren et al. 1996). 

An interesting situation occurs when a muon hops in 
an antiferromagnetically ordered lattice. With no hop- 
ping, and if there is one muon site close to each spin in 
the system, there will be a single precession frequency. 
Because the sign of the internal field reverses every time 
the muon hops, the muon will precess one way and then. 
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after hopping, precess back again. As the hop rate in- 
creases the oscillations are progressively destroyed un- 
til eventually no relaxation is possible (Keren et al. 
1993). This effect has been observed experimentally 
in the cuprate compound Cao.8gSro.i4Cu02 (Keren et 
al. 1993). 



7 Muons and superconductivity 

One of the most fruitful areas of recent research with 
muons has been in the area of superconductivity. The 
last couple of decades have witnessed a renaissance in 
this field following the discovery of high-temperature 
superconductors, organic superconductors, borocarbide 
superconductors, and even superconductors based on 
Cgo buckyballs. To understand the usefulness of muons, 
recall that the two important lengthscales in supercon- 
ductors are the penetration depth, A, which controls the 
ability of the superconductor to screen magnetic fields, 
and the coherence length, ^, which controls the length- 
scale over which the order parameter can vary without 
undue energy cost. If the former is sufficiently greater 
than the latter (the condition is that A > the 
material is a type II superconductor which if cooled 
through its transition temperature, Tc, in an applied 
magnetic field remains superconducting everywhere ex- 
cept in the cores of the superconducting vortices which 
usually are arranged in a triangular lattice. Each vor- 
tex is associated with a magnetic flux equal to one flux 
quantum $o = h/2e. The distance between vortices, 
d, is such that the number of vortices per unit area 
2/^3^^ equals the number of flux quanta per unit area 
B/^o- Thus d oc B^^/^. In general the vortex lattice 
will be incommensurate with the crystal lattice and, ex- 
cept at very high magnetic field, the vortex cores will be 
separated by a much larger distance than the unit-cell 
dimensions. Implanted muons will sit at certain crys- 
tallographic sites and thus will randomly sample the 
field distribution of the vortex lattice. 

In the normal state (T > T^) with a transverse field 
B, all muons precess with frequency lo = (Fig- 
ure |l^(a)). In the superconducting state however the 
muons implanted close to the vortex cores experience 
a larger magnetic field than those implanted between 
vortices. Consequently there is a spread in precession 
frequency, resulting in a progressive dephasing of the 
observed precession signal (Figure |l^(b)). The larger 
the penetration depth, the smaller the magnetic field 
variation and the less pronounced the dephasing (com- 
pare Figure |l^(b) and (c)). It turns out that this idea 
can be quantified (see Aegerter and Lee 1997) and that 
the relaxation rate a of the observed precession signal 




Figure 16: The field distribution inside a supercon- 
ductor as a function of position and the corresponding 
muon-spin relaxation function for three cases: (a) the 
normal state, (b) the superconducting state, (c) as (b) 
but with a shorter penetration depth. 



is related to the penetration depth using 

a = 7^(i?(r) - (i?(r))r)y2 « 0.06097^<i>o/A^ (5) 

where B{r) is the field at position r and the averages 
are taken over all positions. Thus the relaxation rate of 
the observed precession signal can be used to directly 
obtain the magnetic penetration depth. An advantage 
is that data are obtained from the bulk of the supercon- 
ductor, in contrast to techniques involving microwaves 
which are only sensitive to effects at the surface. 

This principle has been applied to many different su- 
perconductors to extract both the penetration depth 
and its temperature dependence. This latter quantity 
is of great interest because it is a measure of the tem- 
perature dependence of the order parameter and can 
yield information concerning the symmetry of super- 
conducting gap and hence the symmetry of the pairing 
mechanism. For example this approach has revealed 
unconventional pairing in a sample of the high tem- 
perature superconductor YBa2Cu3 06.95 (Sonier et al. 
1994). It is also possible to extract the vortex-core ra- 
dius from a detailed analysis of the data (Yaouanc et 
al. 1997, Sonier et al. 1997). 

A conventional type II superconductor exhibits 3 
well-defined phases for T < Tc: (1) a Meissner phase 
for B < Bel, (2) a mixed or Shubnikov phase for 
Bel < B < Bc2 (in which the magnetic field enters the 
superconductor in the form of well defined flux lines or 
vortices arranged in a lattice) and (3) the normal metal- 
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Figure 17: (a) The field distribution p{B) in the vortex 
lattice (contours of B shown in (b)). 



lie phase for B > Bc2- In highly anisotropic systems 
the vortex lattice is no longer a system of rigid rods 
but should be considered as a system of flexible inter- 
acting lines. A useful picture is that of a weakly cou- 
pled stack of quasi-two-dimensional (q2D) "pancake" 
vortices, each one confined to a superconducting plane 
(Blatter et al. 1995, Clem 1991). The phase diagram is 
thus substantially altered to take account of field and 
temperature dependent changes in the vortex lattice 
itself. At low T and low B the stacks resemble conven- 
tional vortex lines. Above a characteristic temperature 
Tb, but still below that at which superconductivity is 
destroyed, the vortex lattice is broken up by thermal 
fluctuations (Clem 1991) (this is called vortex lattice 
melting). At low T, but this time increasing U, the 
energetic cost of interlayer deformations of the lattice 
(local tilting of the lines) is progressively outweighed 
by the cost of intralayer deformations within the super- 
conducting plane (shearing). Above a crossover field 
Bcr the vortex lattice enters a more two-dimensional 
regime. Thus in anisotropic systems we may expect 
temperature and field dependent transitions in which 
the vortex lattice is destroyed. When muons are im- 
planted into a superconductor in a field -Bappiicd one 
can directly measure the field distribution p{B) which 
is given by p{B) = {S{B — B{r)))r and is the proba- 
bility that a randomly chosen point in the sample has 
field B' (Brandt 1988, Aegerter and Lee 1997). This 
is shown in Fig. |lj(a) for an ideal vortex line lattice. 
The distribution is highly asymmetric, the high field 
"tail" corresponding to regions of the lattice close to 
the vortex cores (see Fig. [l7|(b)). The maximum of 
the distribution occurs at i?pk, which lies below the 
mean field (B) (see Fig. |l^(a)). Such lineshapes have 
been observed at low temperatures and fields in vari- 
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Figure 18: The Uemura plot, showing the correlation 
between the superconducting transition temperature Tc 
and the Fermi temperature Tp. The exotic supercon- 
ductors fall within a common band. Adapted from Ue- 
mura et al. 1991 and Hillier and Cywinski 1997. 



ous anisotropic superconductors using /iSR, including 
the high temperature superconductors (Lee et al. 1993, 
Aegerter and Lee 1997) and also in organic supercon- 
ductors (Lee et al. 1997). In both case it is found that 
the vortex lattice can be melted with temperature at 
Tb or can cross into a two-dimensional regime at fields 
above Ba - Both transitions can be followed by measur- 
ing the field and temperature dependence of the p{B) 
line shapes. 

The London formula for the zero temperature limit of 
the penetration depth A(0) yields in the clean limit (the 
mean free path much bigger than the coherence length) 
the relation: A(0) oc y^m* /ns{0) where m* is the effec- 
tive mass and UsiO) is the density of superconducting 
electrons. This can be combined with measurements of 
the Sommerfeld constant to yield a value for the Fermi 
temperature Tp. Thus from muon measurements it is 
possible to plot a diagram showing the relationship be- 
tween the Fermi and critical temperatures for a range 
of superconductors (Uemura et al. 1991). This picture 
(see Figure |l^) has come to be called a Uemura plot 
and shows a clear correlation between Tc and Tp for the 
heavy fermion, organic, fuUerene and Chevrel phase su- 
perconductors. The conventional, elemental supercon- 
ductors lie away from this correlation and have values 
of Tc/Tp ~ 10~^ (for the "exotic" superconductors this 
value is one or two orders of magnitude larger). This 
correlation has been interpreted as evidence that the 
exotic superconductors may be close to Bose-Einstein 
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condensation which is expected to occur at a tempera- 
ture ~ Tp (Uemura et al. 1991). Whether or not this 
speculation is correct, it is expected that this remark- 
able correlation should constrain theories to explain the 
superconductivity in these various exotic systems and 
is possibly suggestive of a common mechanism lying 
behind them. 

8 Muons as active probes 

In almost everything we have considered so far we have 
been tacitly assuming that the muon is a passive probe 
and does not disturb its surroundings. If it has been 
sensitive to dynamics, we have believed that the muon 
takes no part in them itself. This in fact is very of- 
ten true. However, there are a number of situations 
in which the muon plays an active role. For example, 
at high temperature in copper the muon diffuses from 
interstitial site to interstitial site. In this case the ma- 
jor component of the observed depolarization is due to 
the muon motion. In semiconductors the muonium can 
undergo charge and spin exchange with conduction elec- 
trons and thereby one measures dynamics with which 
the muon itself is intimately involved. 

An extreme case where the muon plays a strongly ac- 
tive role is found in conducting polymers (Hayes 1995). 
Figure |l^ shows the reaction between muonium and 
trans-polyacetylene (Nagamine et al. 1984) which pro- 
duces a diamagnetic, neutral muon defect and a highly 
mobile unpaired spin. This soliton diffuses up and down 
the chain but cannot cross the muon defect which acts 
as a barrier. Every time the soliton briefly revisits the 
muon, the muon-electron hyperfine coupling is turned 
on and then off, so that successive visits progressively 
relax the muon polarization. Measurement of the field 
dependence of this relaxation yields the spectral den- 
sity function associated with the defect random walk 
and can be used to infer the dimensionality of the soli- 
ton diffusion (Nagamine et al. 1984). This occurs be- 
cause the relaxation rate is connected with the noise 
power, J{uj^), in the fluctuating magnetic field at the 
muon Larmor frequency, oj^, associated with that par- 
ticular magnetic field. Sweeping the magnetic field al- 
lows one to extract the frequency distribution of the 
fluctuations. In other polymers, such as polyaniline, 
the muon-induced defect is a negatively charged po- 
laron. Muons are uniquely sensitive to the motion of 
this defect in undoped materials (Pratt et al. 1997) and 
in contrast to transport studies, which are inevitably 
dominated by the slowest component of the transport 
process, can provide information on the intrinsic trans- 
port processes governing the mobility of an electronic 
excitation along a chain. 
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Figure 19: Muonium interaction with trans- 

polyacetylene to produce a diamagnetic radical and a 
mobile neutral soliton. 



In more conventional materials it is the motion of the 
muon itself which is of special interest. The dynamics 
of light atoms such as hydrogen and muonium, or par- 
ticles such as the proton and muon, are worthwhile to 
study because they can provide a stringent test for theo- 
ries on the quantum motion of defects and interstitials 
(Storchak and Prokof'ev 1998). The smaller mass of 
the muon leads to larger tunnelling matrix elements to 
neighbouring sites and thus enhances the quantum me- 
chanical nature of the motion. Furthermore, because 
muon-muon (or muonium-muonium) interactions can 
be neglected, the intrinsic nature of the dynamics can 
be followed without the complications that can be found 
in studying the corresponding proton or hydrogen case. 

In inorganic materials the muon will usually come 
to rest at an interstitial site. The stability of that 
site will depend on the depth of the potential well. 
It is of interest to discover whether local diffusion is 
possible between interstitial sites. Another process is 
trapping and release from deep potential wells associ- 
ated with imperfections or defects. The muon jump 
rates are found to be about ten times higher than the 
corresponding proton jump rates, consistent with the 
lighter muon mass. Hopping is, as expected, assisted by 
phonons and thus rises with temperature, following an 
approximately activated behaviour [the hop rates are 
proportional to T~^/^ exp(— i^a/ZcsT) where Ea is an 
activation energy (Flynn and Stoneham 1972)]. This 
occurs because the muon is initially 'self-trapped' by 
its own local distortion of the lattice (Figure 1^) and 
a tunnelling transition is only possible if, by the ther- 
mal fluctuations of the lattice which occur because of 
phonons, two neighbouring energy levels coincide [the 
coincidence configuration, see Figure ^(b)]. The muon 
can then tunnel through the barrier and becomes self- 
trapped in the next site. 

Lower temperature produces fewer phonons and 
hence the hop rate falls as the temperature is reduced. 
However as the temperature is lowered further, a pe- 
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Figure 20: Diffusion process of a muon by phonon- 
assisted tunnelling, (a) The muon is stable in an inter- 
stitial site and the local distortion leads to self-trapping 
so that its zero-point energy level lies a little lower than 
the neighbouring site, (b) Thermal fluctuations provide 
the opportunity for a coincidence configuration whereby 
tunnelling is allowed, leading to (c) a new stable con- 
figuration. (Adapted from Cox et al. 1987). 



culiar effect is observed: the muon hop rate falls to a 
minimum and then begins to rise again. In the very low 
temperature regime the phonons appear to be hinder- 
ing hopping rather than helping it, as they do at higher 
temperatures. The reason is that at low temperatures 
coherent tunnelling is possible i.e. the muon is in a band 
state. Phonons now are responsible for inelastic scat- 
tering which destroys the coherence of this delocalized 
state. This coherent effect is known as quantum diffu- 
sion (for a review, see Storchak and Prokof'ev 1998). 

Experimental data for various materials are shown 
in Figure |2| and although there are large differences 
in the size of the hop rate and the detailed form of 
the temperature dependence, all show an increasing 
hop rate at high temperature consistent with activated 
behaviour and a decreasing hop rate at low tempera- 
ture consistent with quantum diffusion. An early the- 
oretical treatment of this latter effect predicted that 
the low temperature hop rate would follow an inverse 
power law T~", where the exponent a was large, typ- 
ically ~ 9 (Kagan and Klinger 1974). Experiments on 
Cu and Al (see Figure ^l|) showed a more modest be- 
haviour with a ~0.6-0.7. However it was shown that 
this could be explained by considering the effect of the 
conduction electrons in a metal, which screen the muon 
and cannot react fast enough to the diffusing particle 
and follow it adiabatically. This produces a net drag 
which reduces the particle hop rate and weakens the 
temperature dependence. A detailed theory of this ef- 
fect produces agreement with experiment (Kondo 1984, 
Yamada 1984, Kagan and Prokof'ev 1986). 

This dominant role of the electrons has been demon- 
strated by ingenious experiments on aluminium (Karls- 
son et al. 1995). In its superconducting state, the pres- 
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Figure 21: Muon hop rates as a function of tempera- 
ture for various materials. Copper (Luke et al. 1991), 
aluminium (Hartmann et al. 1988), gallium arsenide 
(Kadono et al. 1990), potassium chloride (Kiefl et al. 
1989), and solid nitrogen (Storchak et al. 1994). 



ence of the gap in the electronic spectrum effectively 
decouples the electron bath from the muon. At low 
temperature the superconductivity can be removed by 
applying a sufficiently large magnetic field to the sam- 
ple. This dramatically reduces the muon diffusion rate 
because the closing of the gap reconnects the muon dif- 
fusion process to the electron bath, introducing drag. 

In insulators there are no conduction electrons to 
worry about and in this case the hopping particle is 
a neutral muonium atom, not a charged muon. The 
coherent muon hop rate for KCl (Figure pl| ) rises much 
more rapidly with decreasing temperature than for met- 
als, and fits to an exponent a — 3.3 (Kiefi et al. 1989). 
This still does not quite fit with the earlier theory (Ka- 
gan and Klinger 1974) but agrees with more sophis- 
ticated treatments (see Storchak and Prokof'ev 1998) 
which take into account the measured phonon spectrum 
in KCl measured using neutron scattering. 

In semiconductors muonium is also formed (see sec- 
tion 1^) and very similar temperature dependence is 
found (Figure Kadono et al. 1990) with a ~ 3 be- 
low 100 K but the hop rate saturates below ~10 K due 
to the presence of disorder. (Below this temperature 
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the coherent tunneUing is dominated by the disorder, 
rather than the phonons, and is therefore temperature 
independent.) Data are also shown for sohd nitrogen in 
Figure This material shows a very sharp increase in 
hop-rate with decreasing temperature (a ~ 7, Storchak 
et al. 1994); as the temperature is lowered further it 
saturates and thereafter begins to decrease. This low 
temperature effect is thought to be due to orientational 
ordering of the N2 molecules but is not yet understood 
in detail. 

9 Conclusions 

Muons have found extremely wide application in con- 
densed matter physics, allowing physicists and chemists 
to study a variety of problems concerning magnetism, 
superconductivity, chemical kinetics, diffusion, molec- 
ular dynamics, and semiconductor physics. The inter- 
ested reader is referred to a number of excellent spe- 
cialist reviews for further details concerning the tech- 
nique and its apphcation. (Schenck 1985, Cox 1987, 
Brewer 1994, Karlsson 1995, Dalmas de Reotier and 
Yaouanc 1997). Currently there is also active interest 
in attempting to utilise muons to catalyse fusion (Pono- 
marev 1990). 

The availability of muons at proton sources in vari- 
ous places in the world has led to a great deal of ac- 
tivity in this field. This is because even though the 
proton source is often designed for some other purpose 
(the production of neutrons for scattering experiments 
or the production of mesons for nuclear and particle 
physics), a muon facility can be "tagged on" for little 
additional cost. The consequent reduction in proton 
intensity for the main purpose of the proton source is 
only a few per cent. 

One current limitation of muon techniques is that 
they require bulk samples because the incident muons 
are formed with energy 4 MeV and penetrate a few hun- 
dred micrometres into any sample. Even with the use of 
degraders, surface studies have so far not been possible. 
Recently however there has been considerable progress 
in the development of so-called 'slow muon' beams in 
which the energy of the muon beam is reduced down to 
'^I-IO eV. This can be achieved by either moderation 
in thin layers of rare gas solid (Morenzoni et al. 1994) or 
by resonant ionization of thermal muonium (produced 
from the surface of a hot tungsten foil placed in a pulsed 
proton beam) by a pulsed laser source (Nagamine et al. 
1995). The efficiency of these processes is at present 
rather too low to allow routine experiments on thin 
films to take place. Nevertheless as moderator efficien- 
cies improve and the necessary technology is developed, 
this promises to be a fruitful area of future research. 
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